LA-UR-07-5874 



Calculating the charged particle stopping power 
exactly to leading and next-to-leading order 

Robert L Singleton Jr 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545 USA 
E-mail: bobsl@lanl.gov 

Abstract. I will discuss a new method for calculating transport quantities, such as the 
charged particle stopping power, in a weakly to moderately coupled plasma. This method, 
called dimensional continuation, lies within the framework of convergent kinetic equations, and 
it is powerful enough to allow for systematic perturbative expansions in the plasma coupling 
constant. In particular, it provides an exact evaluation of the stopping power to leading and 
next-to-leading order in the plasma coupling, with the systematic error being of cubic order. 
Consequently, the calculation is near-exact for a weakly coupled plasma, and quite accurate for 
a moderately coupled plasma. The leading order term in this expansion has been known since 
the classic work of Spitzer. In contrast, the next-to-leading order term has been calculated 
only recently by Brown, Preston, and Singleton (BPS), using the aforementioned method, to 
account for all short- and long-distance physics accurate to second order in the plasma coupling, 
including an exact treatment of the quantum-to-classical scattering transition. Under conditions 
relevant for inertial confinement fusion, BPS find the alpha particle range in the DT plasma to 
be about 30% longer than typical model predictions in the literature. Preliminary numerical 
studies suggest that this renders the ignition threshold proportionally higher, thereby having 
potential adverse implications for upcoming high energy density facilities. Since the key ideas 
behind the BPS calculation are possibly unfamiliar to plasma physicists, and the implications 
might be important, I will use this opportunity to explain the method in a pedagogical fashion. 

For a weakly to moderately coupled plasma, the charged particle stopping power dE/dx 
was recently calculated from first principles in Ref. pQ using the method of dimensional 
continuation [2]. While the calculational techniques were imported from quantum field theory 
the calculation itself lies squarely within the standard framework of convergent kinetic equations. 
I will assume some familiarity with Ref. pQ, although a much shorter and self-contained 
exposition can be consulted in Ref. [3j. For ease of presentation, I will work here with a 
one component plasma of charge e, temperature T, and number density n, although the work 
described above is in the context of a general multi-component plasma. Our starting point will 
be the fact that any plasma quantity can be written as a series expansion in integer powers of 
a dimensionless plasma parameter [3]. For example, the stopping power takes the form 

dE e 2 K 

-j— = — A g 2 In g + B g 2 + 0(g i ) , with coupling parameter g = ^° , (1) 

LO NLO error 

where = An e 2 n/T. I have indicated the leading order in g term (LO) and the next-to-leading 
order term (NLO) in Eq. (pQ), while the minus sign in the LO term is a convention that renders A 
positive when energy is transferred to the plasma. Note that nonanalytic terms such as In g can 
also appear in the ^-expansion. For the process of energy exchange via Coulomb interactions, 



this non-analyticity arises from the competition between disparate physical length scales. The 
coupling g is just the ratio of the Coulomb energy, for two particles a Debye length A D = k d ! 
apart, to the thermal kinetic energy as measured by T. To get a feel for the numbers, one finds 
g = 0.04 at the center of the sun, while in the ICF ignition regime g can be smaller. In such 
cases, provided we can calculate A and B, expression (pQ) therefore gives an accurate evaluation 
of dE/dx. 

The coefficient A was first calculated by Spitzer some time ago, while B was recently 
calculated in Ref. pQ using a regularization technique from quantum field theory called 
dimensional continuation. It is convenient to define the dimensionless coefficient C by 
B = —A InC, along with K = Ag 2 , and to then express the stopping power as 

^ = ETlnA coul + 0(g 3 ) , with In A coul = - ln{C g} . (2) 

We see, then, that knowing the next-to-leading order term is equivalent to knowing the exact 
coefficient C under the logarithm. In the extreme classical and quantum limits, the logarithm 
can be written as the ratio of two length scales 

lnA coul = ln{^) , (3) 

where 6 max is set by the Debye screening length k" 1 , and the scale 6 m i n is set either by the 
distance of closest approach in the extreme classical limit, or the thermal De Broglie length 
in the extreme quantum regime. In either case, we find & m in/frmax °c g, and we see that the 
^-dependence in the Coulomb logarithm arises quite naturally. Reference |T] can therefore be 
thought of as a calculation of the Coulomb logarithm, including the exact interpolation between 
the extreme classical and quantum limits. 

Let us now turn to convergent kinetic equations. As suggested by Ref. [5J, one can view 
the Boltzmann and Leonard-Balescu equations as providing complementary physics since they 
both succeed and fail in complementary regimes. The Boltzmann equation (BE) gets the short- 
distance physics correct, while the Leonard-Balescu equation (LBE) captures the long-distance 
physics; conversely, the BE and the LBE miss the long- and short-distance physics, respectively. 
This complementarity motivates a class of kinetic equations of the form [6] 

^ + v . V f = B[f]+L[f]-R[f] , (4) 

where R[f] is a carefully chosen "regulating kernel" designed to "subtract" the long-distance 
divergence from the scattering kernel B[f] of the BE and the short-distance divergence from 
L[f] of the LBE. Furthermore, the kernel R[f] must also preserve the correct physics in the 
complementary regimes, namely, it must not damage the correct short-distance physics of the 
BE and the correct long-distance physics of the LBE. Reference [1] can be viewed as a systematic 
and rigorous implementation of this procedure, albeit in a more abstract form, accurate to second 
order in g (with more work, one could systematically calculate to third and higher order in g). 

We now examine how dimensional continuation regulates the kinetic equations. For simplicity, 
we concentrate on the classical regime, although to second order in g, quantum mechanics can be 
included by using the quantum scattering amplitude in B[f]. The classical BBGKY hierarchy for 
the Coulomb potential is well defined and finite. We run into divergences only when truncating 
the hierarchy to derive lower-order kinetic equations, such as the Boltzmann and the Lenard- 
Balescu equations. Interestingly, this truncation problem occurs for the Coulomb potential, and 
only then in three spatial dimensions v = 3. Therefore, we can regulate the theory, rendering 
it completely finite and well defined, by performing the integrals in an arbitrary number of 



dimensions v. Upon using this procedure, logarithmic divergences in three dimensions become 
finite simple poles 1/(V — 3) in arbitrary dimensions. One can then work entirely with finite 
quantities. In the case of the stopping power, we find that the long-distance pole from the BE 
exactly cancels the short-distance pole from the LBE, and the result is therefore finite when we 
set v = 3 at the end of the calculation. This provides a finite and well-defined result, obtained 
from a regularization prescription constructed in a consistent fashion at all length and energy 
scales. This is a common and time honored regularization procedure in quantum field theory, 
where it is called dimensional regularization. 

I will now review some of the more salient features of the method. Let x„ and v u denote 
the i/-dimensional position and velocity vectors of a particle. The Coulomb potential for two 
particles separated a distance r = |x„ — x'J is V v (r) = C v e 2 /r u ~ 2 , where C v = Y(y jl — l)/4-W 2 
is a spatially dependent geometric factorial The distribution function f v will be defined so that 
fv(x u ,v u , t) d v x d v v gives the number of particles in a small hyper- volume d v x about x„ and d v v 
about Vj, at time t. We can define multi-point correlation functions in a similar manner, and 
in this way we can construct the BBGKY hierarchy in an arbitrary number of dimensions. In 
dimensions v > 3, the standard textbook derivation of the BE goes through without an infrared 
divergent scattering kernel B u [f]. Furthermore, since the i/-dimensional Coulomb potential 
V v {r) oc l/r u ~ 2 emphasizes short-distance over long-distance physics when v > 3, the BBGKY 
hierarchy reduces to the Boltzmann equation to leading order in g in these dimensions: 

df 

BBGKY -j£ + v„- V x f u = B u [f] to LO in g for v > 3 . (5) 

Here, the ^-dimensional spatial gradient has been denoted by V 2 . Conversely, in dimensions 
v < 3, the Coulomb potential V u (r) emphasizes long-distance physics over short-distance effects, 
and consequently, to leading order in g, the BBGKY hierarchy reduces to the Lenard-Balescu 
equation in this spacial regime: 

df 

BBGKY + v„- V x /„ = L„[f] to LO in g for v < 3 , (6) 

where the scattering kernel in the LBE is L v [f}. Space does not permit us to write down the 
exact forms of B u [f] and L u [f] here, but one may consult Ref. [TJ for the expressions. These 
kinetic equations allow one to calculate the stopping power in v > 3 and v < 3, the results 
of which are presented in Sections 8 and 7 of Ref. pQ, respectively. The calculation involves 
performing a series of momentum and wave number integrals in arbitrary dimensions u, and 
reduces to the form 

dF > a 2 

= H(u) + 0{v - 3) : LO in g when v > 3 , (7) 

dx v — 3 

= Q{y) | + 0(3 -v) : LO in g when v < 3 . (8) 

dx 3 — v 

The analytic expressions for H(y) and G(y) are rather complicated!! and space does not permit 
their reproduction here. In this paper, we are only interested in their analytic properties as 
a function of v. In particular, the coefficients H(u) and G{y) can be expanded in powers of 
e = v — 3, and we find 

H{y) = -A + eHi + Oie 2 ) and G{y) = -A + e G x + 0(e 2 ) . (9) 

1 See Ref. 0] for more details. 

2 For the related process of electron-ion temperature equilibration, in contrast, the expressions for H(v) and 
G(y) are quite simple. 



For our purposes, we do not require the exact forms of H\, G\, nor that of the leading term A. 
It is sufficient to note that the leading terms in Eq. ([9j) are equal, so that H{y = 3) = G(y = 3). 
This is a fact that arises from the calculation itself, as it must, and it should be emphasized that 
this equality is not arbitrarily imposed by hand. It is a crucial point that the leading terms are 
identical, as this will allow the short- and long-distance poles to cancel, thereby giving a finite 
result. 

Since the rates dE > jdx of Eq. ([7]) and dE K jdx of Eq. ([5]) were calculated in mutually exclusive 
dimensional regimes, one might think that they cannot be compared. However, even though 
Eq. © was originally calculated in v < 3 for integer values of v, we can analytically continual 
the quantity dE < jdx (viewed as a function of dimension v) to real values of v with v > 3. 
We can then directly compare Eqs. ([7]) and ©. Upon writing the ^-dependence of Eq. ([H]) as 
g2+(u-3) ^ -yyhgn ;/ > 3 we see that Eq. (jSJ) is indeed higher order in g than Eq. (J7J): 

dE< 2+(^-3) 

— — = -G(y) h Oiy - 3) : NLO in g when v > 3 . (10) 

By power counting arguments, no powers of g between g 2 and g v ~ l can occur in Eq. d7J) for 
v > 3, and therefore Eq. (JSJ) indeed provides the correct next-to-leading order term in g when 
the dimension is analytically continued to v > 3, The individual pole-terms in Eqs. dTJ) and (|10p 
will cancel giving a finite result when the leading and next-to-leading order terms are added. 
The resulting finite quantity will therefore be accurate to leading and next-to-leading order in 
g as the v — > 3 limit is taken: 



dE 

—— = lim 

dx 1/^3+ 



dE> dE 
+ 




dx 

NLO 



+ 0(g 3 ). (11) 



Note that this does not lead to any form of "double counting" since we are merely adding the 
next-to- leading order term (|10p to the leading order term (|7|) at a common value of v > 3. We 
are now in a position to evaluate the limit in Eq. (jllj) . Defining e = v — 3 as before, note that 
g e = expjelng} = 1 + elng + 0(e 2 ), which gives the relation 

9 - = - + \ng + 0(e) . (12) 
e e 

Substituting Eq. (fT2"j) into Eq. (fTU|) . adding this result to Eq. dTJ, and then taking the limit gives 

dE 

— = -Ag 2 lng + Bg 2 + 0(g 3 ) , (13) 

with B = Hi — G%, in agreement with Eq. ([1]). In this way, BPS has calculated the charged 
particle stopping power accurate to leading order and next-to-leading order in g. 
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A In the same way that the factorial function n! on the positive integers can be generalized to the Gamma 
function T(z) over the complex plane, including both the positive and negative real axes. 



